Rung— singlet charge— ordering in cV-NaX^Os 
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We show that the proposed zig-zag charge-ordering model for a'-NaV20s is incompatible with a 
single transition that induces both charge-ordering and spin gap. We introduce a two-band model 
for a'-NaV205 within which the simultaneous charge-spin order transition is driven principally by 
electron-phonon interactions. The spin gap is due to the formation of local dimers, thus explaining 
the weak magnetic effect on this transition. 

PACS numbers: 71.30.+h, 71.45.Lr, 75.10-b 



a'-NaV205 has attracted considerable attention be- 
cause of a peculiar insulator-insulator transition at T c 
= 34 K with the simultaneous appearance of charge- 
ordering (CO) and spin gap (SG)1. The material con- 
sists of V-0 ladders coupled through direct V-V bonds, 
with the V-0 layers separated by layers of Na + ions. 
The average number of d-electrons per V-ion is 1/2. The 
insulating behavior above T c is understood if each V-O- 
V rung containing one d-electron is considered a "site" 
within a Mott-Hubbard picture. The SG transition to a 
nonmagnetic insulator 1 was originally explained within 
a "chain model"—, which assumes CO into alternate V 4+ 
and V 5+ chains at T > T c , and spin-Peierls (SP) transi- 
tion in the V 4+ chains at T < T c . NMR demonstrations 
of only one kind of V-ions 3 - at T > T c precludes this 
model. Inequivalent V-ions are, however, found for T < 
T c . 

Subsequently, there have been a number of theoretical 
discussions of the low T state of a'-NaV 2 Oj 4 i 5 i 6 ' 7 ' 8 ' 9 ' 10 , 
also based on the one-band V-only picture. Several of 
these recent theories predict the so-called "zig-zag" CO 
(see Fig.Q)d£, which is driven by the Coulomb repulsion 
between electrons occupying neighboring V ions within 
an extended Hubbard model, 



H = - t rj {dld ja +d) a d la ) (1) 

<ij>,<? 



<y> 



d- creates an electron with spin er on a V 
= d\„di„, riA=nf++nf,, and < ij > implies 



In the above 

d-orbital, nf a = d^cka, rq=n\ 

nearest-neighbor V-ions. The hopping integrals £jj will 
be denoted by t±, t\\ and t a b for rung, leg and interladder 
V-V bonds, respectively. For Uq, U\ » the ground 
state is the zig-zag CO of Fig.^ The origin of the SG is 
not clear, and different theories have proposed different 
bond distortions that generate the SG^&S. 

Experimentally, X-ray diffraction studies agree that 
below T c there occurs a superlattice with 2a x 26 x 4c 
supercell. The period doublings along a and b directions 
and the quadruplings along the c direction must originate 



from lattice displacements. Anomalous X-ray scattering 
experiments have confirmed CO, and have also claimed 
to have detected the zig-zag C O 11 : 12 , although other X- 
ray analysis have suggested that the ordering occurs only 
on alternate ladders 1 ^. Importantly, the relationship be- 
tween the SG and the CO is not obvious from these stud- 
ies, and there is no consensus on the pattern of bond 
distortions leading to the SG. Many experimental inves- 
tigations have focused on proving that the zig-zag CO 
is more appropriate than the chain model, an approach 
that misses the possibility of a third different CO. Op- 
tical measurements indicate a strong role of phonons in 
the CO-SG transition 14 . Magnetic measurements have 
shown that (a) the spin gap 2A S is nearly twice that pre- 
dicted from the BCS relation 2A S = 3.5T,J£, and (b) the 
effect of magnetic fields on T c is much weaker than that 
expected for a simple SP transition 16 . 

There are two objectives of the present Letter. First, 
combining recent theoretical results for 1 / 4— filled band 



system: 



,17.18 



and exact diagonalization studies we have 



performed for finite two-dimensional (2D) clusters we 
show that Eq. (1) necessarily predicts distinct CO and SG 
transitions. Second, we present a two-band model for a'— 
NaV 2 Os within which a single co-operative CO-SG tran- 
sition is explained naturally. Specifically, we show that 
within the two-band model incorporating both electron- 
electron (e-e) and electron-phonon (e-ph) interactions, 
V-O-V ladder rungs are alternately electron-rich and 
electron-poor (see Fig. |2J. The SG is a direct conse- 
quence of the CO, and is due to the formation of local 
singlets on the electron-rich rung V-O-V bonds. We 
point out that there is no contradiction between this pic- 
ture and the Mott-Hubbard description of the T > T c 
state. Local dimers, emphasized in several recent theories 
of correlated electrons 1 ^* 1 ^, were also discussed in earlier 
bipolaron theories of Ti and V-oxides^S and provide a 
new paradigm for SG in 2D systems. 

The fundamental problem with Eq. (1) is that CO 
driven by the spin-independent interaction U± should 
have occurred at higher T where the free energy is dom- 
inated by high-spin excitations, and the SG transition 
should have occurred within the CO phase at lower T. 
This has also been recognized by other investigator a 9 ' 10 . 
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FIG. 1: Schematics of the proposed zig-zag CO. Black and 
white circles correspond to V 4+ and V 5+ , respectively, (a) 
Two ladders, periodic boundary conditions along a and b. 
(b) Single ladder with two neighboring chains, open boundary 
conditions along a. 



Recent numerical work 18 for the single ladder has shown 
that for the zig-zag CO to occur, U\ must be greater than 
a critical value U\ c > 2\t\\\, as in one dimension (ID). 
The magnitude of U\ c increases as Uq decreases from 
infinitjji&, exactly as in ID 1 —. This result, taken together 
with the known result that the ferromagnetic state with 
total spin S = S m ax is described by the Uq — > oo limit 
of Eq. (1), implies that U lc (S = S max ) < U lc (S = 0). 

Hence if U\ > U\ c in a 1 NaV205, CO should occur 

at temperatures comparable to the energy of the ferro- 
magnetic state, while the SG would occur only at lower 
T where singlet formation takes place. 

An alternate scenario, not considered previously, is 
that U\ < U\ c in a'-NaV20s and the CO-SG transition 
is driven by a co-operative lattice distortion. The SP 
transition in ID 1 /4— filled band systems with U\ < U\ c , 
for example, is accompanied by a CO, and has been re- 
ferred to as a bond-charge density wavei 7 -. There are 
two ways lattice-driven SG transition and zig-zag CO 
can occur in a'-NaV205: (i) SP transitions involving the 
zig-zag bonds of single ladders, or (ii) bond distortions 
that involve the entire multiple-ladder system. We post- 
pone the discussion of (i) until after we have presented 
the two-band model, when we discuss why this transition 
is not possible. The co-operative nature of the zig-zag 
CO and bond distortions in the multiple-ladder systems 
has been discussed before^, albeit from a different per- 
spective. First, the difference in the site occupancies can 
drive period 4 bond distortions among the interladder 
bonds 4 '». Next, consider nearest neighbor electron hops 
along the legs of the ladders in Fig.^ In addition to the 
direct hops, there occur "indirect hops" between them 
that involve a V-ion of the neighboring ladder. Since the 
indirect hops involve either a charge-rich ion or a charge- 
poor ion of the neighboring ladder, consecutive leg bonds 
become inequivalent upon CO formation. Similarly, in- 
equivalent leg bonds can drive a dimerization of the in- 
traladder zig-zag bonds in the lattice of Fig. ^b)&. Note 



TABLE I: Exact bond orders for the 16-site lattices in Fig. 
All entries are for lattice of Fig. Q^a) except the zig-zag bonds 
labeled with *, which are for lattice of Fig.^b). 





s z 


=0 


s z 


=4 


Bond 


tab = -05eV 
An=0.238 


tab = -leV 
An=0.191 


tat = -05eV 
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0.0280 
0.0275 
0.0261 
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0.1425 
0.1244 
0.1187 
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0.0169 
0.0472 
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0.0934 
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0.1731 
0.1546 


0.1916 
0.0251 


0.0917 
0.0849 


0.0921 
0.0632 



that the coupled CO-bond distortion necessarily implies 
that these arguments can be reversed, and bond distor- 
tions can be considered to drive the CO^. 

The tendency to bond distortions is measured within 
electronic models from calculations of the bond orders^ 7 -, 
defined as, 

Bi d =J2( d tdja+dld i(T ). (2) 

a 

Differences in the bond orders corresponding to equiv- 
alent bonds indicate spontaneous distortions that would 
occur for + e-ph coupling 17 . In Fig. show the fi- 
nite lattices for which we have done exact diagonalization 
studies. CO is achieved in our finite cluster calculations 
by adding a site—energy term £\ efnf to Eq. (1), where 
ef is negative (positive) for the charge-rich (charge-poor) 
V-ions of Fig.^ and can originate from e-ph coupling 9 . 
Our calculations are fo:»2I JJq = 1 — 4 eV, t± = 0.35 eV, 
t|| = 0.15 eV and t ab = 0.05 - 0.2 eV, and < U x < 0.5 
eV. Since the CO and bond order differences are driven 
by the site energies in our calculations, our results are 
qualitatively the same for both Ui = and U\ ^ 0, and 
we report the results for U\ = only. In Table [I] we 
present the results for Uq = 3 eV and \ef\ = 0.1 eV for 
two values of t a h. For both values of t a h there occur pe- 
riod 4 bond distortions along the interladder bonds and 
period 2 distortions of the ladder leg bonds in the S = 
state. In addition there occur period 2 bond distortions 
along the zig-zag bonds for the lattice of Fig. ^b). As 
mentioned above and discussed extensively in reference 
ITtI it is not necessary that the CO drives the lattice dis- 
tortion. Rather, the bond distortions of Table |U can also 
drive a zig-zag CO. 

To investigate one versus two transitions, we have stud- 
ied total spin S = S max = 4 and these results are also 
reported in Table [I] It is seen that the tendency to bond 
distortion is probably even stronger in S = S ma x than in S 
= 0. Note in particular that the difference in the charge 
densities between the charge-rich and the charge-poor 
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FIG. 2: Two-band models for (a) MX chains, (b) single ladder 
of a'-NaV20s. In (b) the large (small) circles correspond 
to 0(V)-ions. Occupancy by electrons is larger (smaller) in 
the black (gray) smaller circles in both (a) and (b). Singlet 
formation occurs on the electron-rich rungs in (b). 

sites, An, is larger in the S = S max state. We have per- 
formed similar calculations also for the case of CO on 
alternate ladders^. Once again, An and the bond dis- 
tortions are larger in the ferromagnetic state than in the 
singlet state. There is therefore no reason for the coupled 
CO-bond distortion to occur only in the S = state. We 
therefore conclude that even with nonzero e-ph coupling 
the zig-zag CO does not explain the observed CO-SG 
transition in a'-NaV20s. This is understandable: ex- 
actly as there are two chains per ladder, there are two 
possible zig-zag CO patterns, and there is no fundamental 
difference between the chain and zig-zag models, except 
that a short range U\ instead of long range Madelung 
interaction has been invoked to stabilize one zig-zag CO 
over the othe»^. 

We now show that a different CO pattern emerges 
when we include the O 2- ions. The most natural two- 
band model for a'-NaV20s is the 2D extension of the 
Hamiltonian that has been successfully used to describe 
ID transition metal-halogen (M-X) chains 2 ^, 

H = Hi,2 + «A JM , v ](4^^ + ^c-) (3) 

<ij>,fi,a 

IS 

<jij'>,H,v <ij>,fi 
<fj.fj.'>,<ij>,a i 

In the above, i, j can be both V and O; d\ and p\^ 
create electrons in V d-orbitals and O p-orbitals, respec- 
tively; fj, is the ladder index; (jij') implies consecutive 
sites in any direction, and {/ifi') implies neighboring lad- 
ders. The V-0 hopping integrals t\ and t^ are along the 
ladder legs and rungs, respectively, and ti are the site 
energies (ey = +e, eo = — e); Ui^ are the displacements 
of the V or O ions from their equilibrium positions and 
Aj Mi i M = — Uj^, a and (3 are the intersite and on-site 
e-ph coupling constants, and K is the spring constant. 
For simplicity we have chosen a, p and K to be same 
for leg and rung bonds. The number operator cor- 
responds to both V and O. We have included the on- 




FIG. 3: The amplitude of the ladder- leg O-ion distortion 
versus the Hubbard repulsion for V-ions for electron-phonon 
coupling constants X a = a 2 /Kti = 0.44, \p = /3 2 jKt\ = 0.11 
(solid curve) and A a = 0.28 \(j = 0.07 (dashed curve). Inset 
shows the behavior for Uo = 1 eV 



site Coulomb interaction only for V. Because of the large 
negative site energy for O, the Hubbard interaction for 
O-ions makes little difference in the final results. We 
have not explicitly included intersite Coulomb interac- 
tions, the effects of which we discuss later. 

We have calculated self-consistent solutions to the Uq 
= limit of Eq. (2) for both single ladders and cou- 
pled multiple-ladder systems, with periodic boundary 
conditions along a and b directions for systems up to 
320 atoms. Typical parameters were t\ — 0.9 eV, t^ — 
1.3 eV, and ey (eo) = 2.5 eV (- 2.5 eV), which were 
derived from reference |^. Stable solutions for (3 = + 
correspond to breathing mode distortions of the O-ions 
shown in Fig. 2(b) for the single ladder, leading to (i) 
inequivalent charges on the V-ions on consecutive ladder 
rungs, (ii) shorter V-O-V bonds on the rungs contain- 
ing nominally V 4+ ions, and (iii) period doubling along 
the b direction. As there are two possible phases be- 
tween adjacent ladders, period doubling in the a direc- 
tion is possible as well 9 . As indicated in Fig. 2(a), there 
is a one-to-one correspondence between the distortions 
in the V-0 ladders and the M-X chains^. Each lad- 
der rung is analogous to a single M-site, with internal 
structure that leads to on-rung bond distortion. There 
also emerges a strong reason to ignore the SP transition 
involving the zig-zag bonds within the single ladder (see 
above). The distortion in Fig. 2(b) corresponds to al- 
ternating leg 0-0 distances and inequivalent V charges; 
the SP distortion of the intraladder zig-zag V-bonds re- 
quires alternating leg V-V distances and nominally O 1- 
and 3 ~ anions. The impossibility of obtaining 3 ~ pre- 
vents this particular distortion. This result is identical to 
the observation that M-M distances are uniform within 
the two-band model for M-X chains 23 . 

The above two-band Uq = calculation is signifi- 
cant, as the shorter bonds on the rungs containing the 
charge-rich V-ions suggest singlet formation. CO and 
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SG therefore occur simultaneously within this model. 
The Uq = model is, however, unsuitable for the in- 
sulating state observed at T > T c . A correct description 
of a'-NaV205 then requires nonzero but small Uq within 
Eq. (3). Specifically, Uq should be large enough that 
the Mott-Hubbard picture for insulating behavior for T 
> T c is qualitatively valid, and small enough that the 
rung-based CO is not precluded. Note that the repul- 
sion between two electrons on the same rung bond is an 
effective Hubbard interaction U e ff that is much smaller 
than the bare Uq, while for small t a h, the effectively ID 
1/2-filled band ladder is insulating for all U e ff. For our 
model to be valid there must exist a range of Uq over 
which the distortion of Fig. [3 is nonvanishing and its 
amplitude decreases slowly. 

We have calculated the O-ion distortion \uo\ for 
nonzero Uq within Eq. (3) for the single ladder with 4 
rungs, using the constrained path quantum Monte Carlo 
technique 2 -. We retain the leg O-ions, but consider only 
effective V-V bonds on the rungs, as the rung oxygens 
play a weak role in the CO and influence only the in- 
trarung distortion. For different a, (3 and K, we calculate 
the distortion \uo\ at which the energy minimum occurs 
for each Uq (see insert, Fig. |3J). The overall results for 
two different sets of e-ph coupling constants are shown in 
Fig. 01 As expected, \uo\ depends strongly on a and /?, 
and also decreases with Uo- Significantly though, while 
the decrease in \uo\ is moderately rapid at small Uq, the 
decrease is rather slow for Uq > 2 eV. This behavior is 
a signature of the persistence of the rung-based CO for 
realistic Uq ~ 3 — 4 eV. 

A strong test of the rung-based CO model involves 
demonstrating the co-operative nature of the CO-SG 
transition. The calculated An and bond distortions in 
Table U for the zig-zag CO are larger for S = S max than 
for S = 0. We have done the same one-band calculation 
also for the rung distortion, for the lattice of Fig. 1(b), 
and for isolated 6 and 8 rung ladders. The model Hamil- 
tonian is again the Hubbard Hamiltonian (Eq. 1 with Ui 
= 0) including the site-energy term ^ efnf , but now the 
ef are equal for V-ions on the same rung, and alternate 
in sign for consecutive rungs. The results are identical in 
all cases, viz., (i) bond orders are larger on rungs with 
larger charge density in the S = state, and (ii) the CO 
and bond distortion is absent in S = S max . For the same 
Uq and ef as in Tabled for example, An ~ 0.1 and ABij 
~ 0.05 in the S = state of the single ladder with 8 
rungs, while both An and AB^ are ~ in the ferro- 
magnetic state. In contrast to the zig-zag model there- 
fore, there occurs a single insulator-insulator transition 



involving both charge and spin within the rung distortion 
picture. While the zig-zag CO could be obtained in a 
two-band model with oxygen breathing distortions out 
of phase on the ladder legs£, this would still give distinct 
CO and SG transitions. 

We now discuss the role of intersite Coulomb inter- 
actions, which were ignored in Eq. 3. Intersite V-0 
Coulomb interactions have a weak effect on the CO-SG 
transition, as might be expected. We have already ar- 
gued that the V-V Coulomb interaction is smaller than 
the critical value necessary to promote the zig-zag CO. 
On the other hand, the second neighbor V-V interac- 
tion, ignored in Eq. 1, strongly destabilizes the zig-zag 
CO and promotes the rung CO. There exists therefore a 
wide range of parameters for which the rung-based CO 
dominates over the zig-zag CO. 

In summary, we have shown that the zig-zag CO is 
not compatible with a single CO-SG transition in a'— 
NaV205, and have proposed a two-band model for this 
system within which a co-operative transition involving 
both charge and spin is driven by breathing mode vi- 
brations. Strong d-p hybridization exceeding the energy 
scale of Uq in V-oxides has been suggested by Zimmer- 
mann et al., based on photoemission studies 2 ^. Interest- 
ingly, the bond-order structure proposed in ref. (see 
Fig. 4) is consistent with rung ordering. Our theory 
is consistent with the observed strong phonon effects^, 
while local tightly bound singlet dimers would be in 
agreement with the observed weak effect of magnetic field 
on the CO-SG transition 16 . Careful reexaminations of 
the X-ray data are called for. As mentioned already, 
many of these experiments were carried out to distin- 
guish between the chain and zig-zag CO's, which are the 
only CO patterns possible within the one-band model. 
Our results demonstrate that the tendency to form lo- 
cal dimers is particularly strong at or near 1/4-filling 
of the bancUiSi. Local dimer formation provides a new 
mechanism for SG transitions beyond the more common 
mechanisms involving SP dimerizations in ID and spin 
frustrations, and it has even been suggested that local 
dimers may be relevant in the context of superconduc- 
tivity in correlated electron systemsi2ii^. It is highly 
interesting in this context that (a) in superconducting 
/3-Nao.33V205 there occur crystallographically three dif- 
ferent kinds of V-ions, with perhaps one class of V-band 
1/4-filled 2 ^, and (b) the same bandfilling characterizes 
organic superconductors. 
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